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ABSTRACT 

A major hindrance to the use of algebraic techniques for the 
design of linear control systems Is the non-llnearlty of the equations 
describing the system* The equations for a general class of control 
systems are shown to be linear In certain groupings of variables 
contained In, or added to, the systems. A design technique based 
on the solution of these linear equations Is outlined, and a 
digital computer program which will accomplish the design Is pre- 



sented 



DUDLEY KWQX LIBRARY 

NAVAL POSTGRADUATE SCHOOL 

MONTEREY CA 93943-5101 Table of Contents 

Section Page 

1. Introduction 11 

2. Algebraic Methods 13 

Method of Mltrovlc' 14 

Coefficient Matching 15 

3. Algebraic Solution of a General Problem 18 

Steady-State Error Specification 24 

Bandwidth Specification 25 

Multiple Specifications 26 

4. Limitations of the Method 29 

Steady-State Error for a Ramp Input 29 

Multi-path Feedback 30 

5. Sample Solutions 32 

6. Conclusions 75 

APPENDIX Program COMPTECH 77 



3 



List of Tables 



Table 

1 . 



Page 

Matrix Systeta of Equations for 

Compensated Feedback Control System. 23 



5 



List of Illustrations 



Figure 

1, C and 0 )^ in the s -Plane 

2. Basic Feedback Control System 

3* Basic Feedback Control System with 
Cascade Compensation 

A. Control System with Two Feedback Paths 

5. Block Diagram for Example 5.1 

6. Block Diagram for Example 5.2 

7. Block Diagram for Example 5.3 

8. Block Diagram for Example 5. A 



Page 

13 

18 

20 

30 

32 

A1 

52 

63 



7 




%• 




¥ 




X(8) 

Xn(s) 

Xd(s) 

XOk 

Kt 

T(s) 

H(s) 

C(s) 

E88 

^db 

U) 

C 

lii 

n 



Table of Symbols and Abbreviations 
a complex variable (8 =» a + jt*!) 

a transfer function; a ratio of polynomials In the 

variable s. 

the numerator of X(s) 

the denominator of X(s) 

the coefficient of the kth power of s in Xn(s) 

a function of the variables x, y and z 

variable forward gain 

forward transfer function 

variable feedback gain 

feedback transfer function 

filter compensation transfer function 

steady>state error 

gain in decibels specified in a bandwidth requirement 
frequency specified in a bandwidth requirement 
damping ratio, cosine of an angle in the s»plane 
undamped natural frequency, a radial distance in the 
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1. Introduction 



Linear control system design may be generally described as the 
synthesis of physical components or systems, active or passive, whose 
dynamic response will satisfy specified criteria. Present design 
procedure Is essentially mathematical, since the response of dynamic 
systems Is conveniently described by a set of differential equations. 
The problem Is then to: 

1) formulate differential equations whose solution will 
describe the desired response, 

2) relate coefficients of the equations to the parameters of 
a physical system, 

3) obtain a solution of the differential equations In terms 
of these parameters. 

Production of systems designed In this manner Is frequently economi- 
cally unfeasible, since components required may be expensive to man- 
ufacture In limited numbers. A more practical approach Is to utilize 
readily available components, combining them and adjusting their 
variable parameters to obtain a system with the desired dynamic 
response. 

Analysis of systems thus synthesized may be accomplished In the 
time domain by direct solution of the system differential equations, 
but the more usual approach Is to perform analysis In the complex 
frequency domain (s-plane), by means of the Laplace Transformation. 
System equations thus transformed may be approached with several well- 
known analysis techniques, notably the Evan's root locus procedure, 
or In the frequency domain (s “ Bode, Nyqulst or Nichols 

methods. Application of these techniques to the design procedure Is 
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restricted, however, since they are not easily described analytically, 
they are essentially graphical and therefore time-consuming, and the 
approach to design is largely trial and error. All suffer the further 
disadvantage that changes in system response may be predicted on the 
basis of variation of only one parameter. 

Several recent techniques, to be discussed in further detail, 
have succeeded in overcoming the more obvious difficulties listed 
above. The transformed differential equations describing the system 
and the specifications prescribing system response are reduced to a 
set of algebraic functions of the variable parameters. The algebraic 
relationships produced, however, are non-linear, and for complex 
systems may prove difficult or impossible to solve. This paper out- 
lines a numerical procedure for the solution of these non-linear 
equations, and presents a digital computer program which obtains 
solutions for a general class of control device problems. 
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2 . Algebraic Methods 

Application of algebraic techniques requires the reduction of 
the transformed system equations to a "transfer function," usually 
the ratio of two polynomials in the variable s, which relates the 
output of the system to an arbitrary input. The roots of the denom- 
inator polynomial, known as the "characteristic equation," determine 
the stability and response time of the system, factors of principal 
concern to the designer. For the system to be stable, i.e. finite 
response to a finite input, it is necessary that all roots of the 
characteristic equation have negative real parts. It can be shown 
that the time response of the system is determined by the values of 
the roots, and is frequently dominated by the root or pair of roots 
whose real parts have the least negative value, i.e. when plotted in 
the s-plane, the roots located closest to the origin. The dominant 
roots are designated in polar coordinates by their radial distance 
from the origin (03^), and the cosine (5) of the angle formed by the 
second quadrant radius vector and the negative real axis. 




5 and w in the s-Plane 
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Method of Mitrovic' 



The characteristic equation is a pol 3 naomial of the form 
F(s) = + F^_^ + + F^s^ + Fj^aV F^ (2.1) 

or in sununation notation 

" k 

F(s) = E F s (2.2) 

k=0 

where the polynomial coefficients (F^^'s) are functions of the system 
parameters. Mitrovic* [1] has shown, through a transform of variable 
and s-plane mapping procedure, that two coefficients of F(s) may be 
stated as functions of the remaining coefficients and of the ^ and 
of a pair of specified roots. The result of this procedure is two 
equations of the form 

E (-1)^ F^a,„^ 4, (O = 0 (2.3) 

k=»0 

? (-1)^ Fj^ <|)^^(0 = 0 (2.4) 

k=0 

where recurring Chebishev functions, which for 

a specified 5 may be precomputed from the relationship 

‘*’k+l^^^ = 0 (2.5) 

with 4> (5) = 0 

o 

and <l>j^(5) = -1 

If numerical values are assigned to all but two of the coeffi- 
cients (F. ) and to 5 and o) , equations (2,3) and (2,4) are reduced 

K. ^ 

to a pair of simultaneous equations in two unknown variables, i,e, 
the coefficients which have not been specified. Mitrovic* originally 
chose the coefficients F^ and Fj^ to be variables, arriving at the 
relations 
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( 2 . 6 ) 



o 

It 


f^ (Fj^, ?, (2.6) 


II 


h <2.7) 



Later, Slkjak generalized the Mitrovic' method to show that a solution 
can be obtained for any pair of coefficients [2] . Thus, a relationship 
between the coefficients of the characteristic equation, which are 
functions of the system parameters, and the location of a desired pair 
of roots is established. Designating the variable system parameters 



as (a^, 


.... a ) , we have for any two coefficients 
z 


F » 

X 


®2* """ ®z^ “ 


II 


*2* ””” * ^6^^k* ^ ’ ^n^ (2.9) 



The process may be repeated to specify additional pairs of roots, 
thus establishing additional relationships of the form of (2.8) and 
(2.9). If, however, not all roots of the characteristic equation 
are specified, the method provides no information on the location of 
undesignated roots, and does not guarantee either stability or the 
dominance of a desired pair of roots. 

Coefficient Matching 

If all of the roots of F(s) are designated, the characteristic 
equation may be stated as 



F(s) « 


n 

® 7 T + r, ) (2.10) 

k«l 



where (r, , r«, r ) are the roots of F(s). The resultant coeffi- 
' 1* 2* n' 

dents of the expanded polynomial are functions of the roots. 



F(S) : 


(2.U) 
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where — - r^) is the sxan of the products of all possible 

combinations of (r, ,r_, r ) taken k at a time. 

1 z n 

For example: 

= (rjrjrj— r^) + (r^rjr^— + (r^Vn-kH-r" V 

+ — Vl> * — - 

In summation notation 
X n 

M®) “ ^ I ^ k (rJl ; X = __nj (2.13) 

i=l j=i ^ kl (n-k) I 

Since all roots are designated^ the Fj^'s are fixed constants. The 
result is two expressions for F(s) 



It 

(0 


i 

I £ (a ,a — a )8 


Parameters 


(2.14) 




i=0 






F(s) = 


" i 

1 g.(r ,r — r )s 

1=0 ^ ^ ^ " 


Roots 


(2.15) 



In order that these expressions be identically equal, coefficients 
of the same powers of s must be equal. 

Fj^s^ = fj^(aj^,a2,— a^)s^ = ®k^^l’^2’ (2.16) 

The resultant system is composed of n simultaneous equations in z un- 
knowns. Complete solution of the system of equations requires that 
the munber of unknown quantities (z), the variable parameters, be equal 
to the number of equations (n), the order of the characteristic equa- 
tion. If insufficient variables are available, additional parameters 
must be inserted in the problem in order that this requirement be 
satisfied. 

The principal advantage of designation of all roots is the guar- 
antee of stability and dominance of desired roots. By designating a 
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pair of roots whose values give the desired time response, and loca- 
ting the remaining roots in such a way as to Insure dominance of the 
original pair and stability of the system, the designer completely 
specifies the time response of the system. The principal disadvant- 
age is the necessity for addition of arbitrary parameters sometimes 
required to permit solution of the simultaneous equations. 

A disadvantage common to all algebraic methods is the Inherent 
non-linearity of the resultant equations. The expressions usvially 
contain products and/or powers of the variable parameters, and 
assistance of a digital computer is required for solution of all 
but the most simple cases. 
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3, Algebraic Solution of a General Problem 



Given a dynamic device whose transformed time domain equation 
may be expressed as the product of a variable gain (Kt) and a trans- 
fer function T(s), and a feedback device expressible as a variable 
feedback gain (Kf) and the transfer function H(s), the transfer 
function F(s) of the "closed loop" control system comprising the two 
components, as arranged in Figure 2, is 



Output (s) 
Input (s) 



F(s) 



Kt T(s) 



1 + Kt T(s) Kf H(s) 



(3.1) 




Figure 2 

Basic Feedback Control System 
Both T(s) and H(s) are the ratio of polsmomials in S. 



Trs'k - 

Td(s) 



Tn s® + Tn + + Tn-S^ + Tn^s + Tn 

a a«*l 2 1 o 

Td.s^ + Td, + + Td-s^ + Td.s + Td 

D b-1 Z 1 o 



^ H- Hn^a H- Hn^) 

Hd s^ + Hd s + Hd 
2 1 o 



(3.2) 



Kt Tn(s) 

and F(s) = Td(s) (3.3) 

1 + Kt Tn(s) Kf Rn(s) 

Td(s) Hd(s) 

F(s) reduces to the ratio of two polynomials 
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F(s) 



g"Cs) 

Fd(s) 



Kt Tn(s)-Hd(s) 

Td(s)*Hd(s) + Kt Kf Tn(s) •Hn(s) 



(3.4) 



Forming the product of polynomials where indicated defines: 

Tn Hd(s) = Tn(s)*Hd(s) 

Td Hd(s) = Td(s)‘Hd(s) 

Tn Hn(s) « Tn(s)»Hn(s) 

The characteristic equation of the system - Fd(s) is the sum of the 
polynomials TdHd(s) and KtKfTnHn(s). 

Fd(s) = TdHd(s) + KtKfTnHn(s) 

= TdHd + TdHd ,s^”^ + + TdHd.s + TdHd (3.5) 

X x-1 1 o 

+ KtKf(TnHnyS^ + TnHny^j^s^”^ + + TnHn^s + TnHn^ 

In most devices used for control systems, particularly higher order 
systems, the order of the numerator polynomial is much less than the 
order of the denominator. Feedback devices seldom have numerator 
order more than two greater than denominator order, i.e. second 
derivative feedback. Thus the order of the characteristic equation, 
and the number of roots to be designated, is determined by the order 
of the TdHd(s) product. By expanding the designated roots into a 
polynomial R(s), and equating coefficients, the system of simultaneous 
equations for Fd(s) becomes 



(WHdjj ) 


X 

8 


= 


R 8 

X 


(WHd^.l) 

• 

» 




S3 


E 

x-1 


♦ 

(TdHd 

y 


+ KtKfTnHn )«^ 

y 


Si 


R sy 

y 


(TdHdy_^ 


+ KtKfTnHny_j^)8^"^ 


Si 


R syi 


• 

(TdHd^ 


+ KtKfTnHn^ )8^ 


= 


R 8^* 
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R.s 

L 



(TdHdj^ ■= KtKfTnHn^ )s 

(TdHd = KtKfTnHn ) R (3.6) 

If all polynomials axe rationalized , 

TdHd = R = 1 

X X 

and the system consists of x equations in the two unknown parameters 
Kf and Kt. In order to obtain an algebraic solution, an additional 
(x-2) variables must be inserted in the system. The addition of 
"cascade compensation" in the form of a multi-section filter net- 
work permits introduction of the necessary mathematical variables 
and is also an economical method of inserting variable parameters 
in the physical plant. 

Cascade Compensation 




Figure 3 

Basic Feedback Control System with Cascade Compensation 

The transformed output/input equations C(s) of a resistance- 

capacitance filter network are expressible as the ratio of two 

equal order polynomials, the order being equal to the number of 

filter sections employed. 

Cn s'" + Cn ,8®"^ + — + Cn,s^ + Cn.s + Cn 
m / I o 

Cd s'" + Cd ,8™“’ + — + Cd. s' + Cd.s + Cd 
m ra-1 ^ 1 o 



C(s) 



Cn(s) 

Cd(s) 



(3.7) 
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The transfer function of the compensated closed loop system, for 
the configuration shown in Figure 3, becomes 



= Y(a\ » Kt T(»)*C(s) 

Input (s) ^ 1 + Kt T(s).C(s) Kf H(s) 

which reduces to 



(3.8) 



F(s) 



Fn(s) _ Kt TnHd(s)>Cn(s) 

Fd(s) “ TdHd(s)*Cd(s) + KtKf TnHn(s)«Cn(s) 



The order of the characteristic equation is now determined by the 
TdHd(s)»Cd(s) product, thus the addition of each filter stage 
Increases by one; the order of Fd(s), the number of simultaneous 
equations to be solved, and the number o£ independent variables 
required. If the order of TdHd(s) is (x), and (m) filter sections 
are used, the order of Fd(s) will be (x + m). The addition of each 
filter section (1) Introduces two Independent variables into the 
system, Cn^ and Cd^. Thus (m) filter sections Introduce (2m) 
variables. 

For solution of the system of equations, the product of the 
Independent parameters (KtKf) may be considered as a single 
algebraic variable. Addition of (m) filter sections to a system 
whose TdHd(s) is of order (x), produces (x + m) equations in 
(2m + 1) variables. Equating the number of variables and number 
of equations 

X + m “ 2m + 1 

m X - 1 (3.10) 

For the system given in Figure 3, the number of filter sections 
required is determined by the product of the denominator polynomials 
of the forward and feedback path transfer functions. 
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The system of simultaneous equations for the compensated system 
is obtained by designating the required (x + m) roots, expanding the 
roots to form the polynomial R(s), forming the TdHd(s)»Cd(s) and 
KtKfTnHn(s)»Cn(s) products, and equating coefficients. Table 1 
Illustrates, in matrix form, the equations for a system having a 
third order Td(s)*Hd(s) product, a first order Tn(s)*Hn(s) product, 
and the required second order filter. The format of the matrix 
equations is the same for higher order systems if the requirements 
that 

1) The number of filter sections used be one less than the 
order of Td(s)*Hd(s), 

and 2) The order of Tn(s)*Hn(s) be less than the order of 
Td(s)‘Hd(s) 
are fulfilled. 

Note that the equations for systems of this type are linear in 
certain individual parameters (Cd^, Cd^^ — - Cd^), and in groupings 
of the other parameters (KtKfCn^, KtKfCn^^, KtKfCn^). Thus the 
system of equations may be solved by linear methods for these 
individual and grouped parameters. 

Again, if all polynomials are considered rationalized. 



TdHd Cd 
X m 



xjfm 



so Cd^ “1. If Cn^ is defined as equal to 1, the group variable 

KtKfCn^ becomes simply the product KtKf. Having solved the linear 

system and obtained the coefficients of the denominator polynomial 

(Cd^, Cdj^, — Cd^_j^) of C(s) directly, the coefficients of the 

numerator polynomial may be obtained by dividing the grouped variables 

(KtKfCn , KtKfCn. , — KtKfCn ,), by the value of KtKfCn , i.e. 
o* V m-1 * ^ m* 
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by KtKf. With the numerator and denominator polynomials 
available, the filter design becomes a problem in passive network 
synthesis. 

A solution has thus been obtained for all physical parameters 
except the individual gains, Kt and Kf, An additional specification 
is required in order to obtain an independent solution for values 
of the gains. The additional criteria considered are steady state 
error and bandwidth. 



Steady - State Error (Ess) 

The error of a system is defined as the difference between 
an input command signal and the system response. 

Error (s) = Input (s) - Output (s) (3.11) 

From Equation 3.8 

Output (s) = Kt T(s)»C(s) _ ^ / N 

1 + Kt T(s)‘C(s) Kf H(s) 

The error as a function of the input signal becomes 

Error (s) - 1 T ( sl »C( sJI [K f H Cs ) , I] . . .3 . 

rror is; 1 + Kt T(s)*C(s) Kf H(s) ^ 

xdiich reduces to 

Error (s) _ TdHd(s)»Cd(s) + KtKf TnHn(s)*Cn( 3 ) - KtTnHd(s)«Cn(s) 



Input (s) 



TdHd(s)#Cd(s) + KtKf TnHn(s).Cn(s) 



(3.13) 



For a unit step displacement input. Input (s) is 1/s. Error (s) 
for a step input becomes 



Error (s) « 






1 - 



Kt TnHd(s)»Gn(s) 



».Cn(s)] 



TdHd(s)»Cd(s) + KtKf TnHn(s) 

Steady-state error is defined as the error of the system after all 
transient effects have taken place, or mathematically, as time 
expired after application of the input signal approaches infinity. 
Application of the Final Value Theorem to Equation 3.14 gives 



14) 
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Ess = 11m 

S-* o 

Ess = lim 
s-*o 



8 Error (s) 



Q 



Kt TnHd(s)»Cn(s) 

TdHd(s).Cd(s) + KtKf TnHn(s)*Cn(s 




TdHd . Cd + TnHn KtKfCn - TnHd . KtCn 
o o o o o o 

TdHd • Cd + TnHn • KtKfCn 



(3.15) 



o o o o 

Note that steady-state error as expressed above Is defined for unity 
Input. Error Is relative to Input, and maximum steady-state error 
desired In the system Is specified as a fraction of the Input signal. 

Steady-state error Is found to be a function of two variables 
previously defined, Cd^ and KtKfCn^, and an additional variable, 
KtCn^. Values for Cd^, KtKfCn^ and KtKf having been previously 

obtained, the solution for Kt Cn Is 

* o 



Kt Cn 

o 



(1 - Ess) [TdHd^»Cd^ + TnHn^. KtKfCn^] 
TnHd 

o 



(3.16) 



and Kf 



KtKfCn 

o 

KtCn 



Kt 



KtKf 

Kf 



(3.17) 



o 

Values of all unknown parameters have been obtained, and a system 
Is synthesized which will meet specified requirements for stability 
root dominance, and steady-state error. 






Bandwidth (BW) 



The bandwidth of a system Is defined as the frequency range 
over which the response of a control system Is within acceptable 
limits [1], It Is frequently given as that frequency (to) at which 
the magnitudes of the output response and the Input signal have a 
specified ratio (G), usually expressed In decibels (db). 

If F(s) Is the closed loop transfer function of a control 
system, bandwidth may be found from 
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20 log 



( output ( .0)) 



10 \ input 



L 



db 

where to is the desired bandwidth 



- 20 log^Q[F(jO))] 



( 3 . 18 ) 



F(j0)) 



Fn(j co) 
Fd(^to) 



Fn ( . oj) + + Fn_(.w) + Fn, + Fn 

a^l 2^i 1 1 o 

Fd^(^w)^ + + Fd2<j(^)^+ Fdj(^w) + Fd^ 



Fn( (i Qreal + j FniLajjlmaf 
Fd(o))real + j Fd(o))iiaag, 

| ~(Fn((i))real)^ + (Fn((jj)jgtag.)^ 

I 2 2 

j_(Fd(a3)real) + (Fd((^)imag.) 



( 3 . 19 ) 



The bandwidth expression is seen to be highly non-linear, however, 
this obstacle is not insurmountable. Recalling that 



F(j«) 



Kt TnHd(^ to)« Cn(^ to) 



TdHd(ja))«Cd(ja)) + TnHn(ja))#KtKfCn(jO)) 



and that values for the variables Cn, , Cd., and KtKfCn, are available 

i i i 

from the solution of the simultaneous equations, Kt may be obtained 
from the bandwidth criterion. If gain (G) is specified at frequency 



(“)» 

Kt 



c 



TdHd(. to)«Cd(^g)) 



TnHd(j ui) »Cn(j oj) 



TnHn( ^ to)* KtKfCn( . 



] 



and Kf - KtKf 
Kt 



( 3 . 20 ) 

( 3 . 21 ) 



Again, values of all unknown parameters are obtained, and the system 
design will meet stability, root dominance and bandwidth specifications. 

Multiple Specifications 

Either a steady-state error or a bandwidth specification is 
sufficient in itself to permit solution for all variable parameters 
and design of the system. However, if both error and bandwidth 
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are specified, solution for one does not guarantee satisfaction of 
the other requirement. A slightly different approach must be taken. 

No requirement has been placed on the location of the non- 
dominant, or non-speclfled roots, other than that they have negative 
real parts, and’ that they be sufficiently removed from the origin of 
the s-plane to insure dominance of the desired roots. This gives the 
designer a great deal of latitude in the choice of extraneous roots 
which are added to the problem as a result of cascade compensation. 

The system of linear equations used in the solution for para- 
meters is obtained by equating coefficients of the characteristic 
equation. Ott one side of each equation is an algebraic function 
of the variable parameters; on the other side is the niimerical 
coefficient found by expanding the roots to form the polynomial 
R(s). Changing the location of one or more roots changes the values 
of the variable parameters. Parameter values obtained are those 
which will place the roots of the characteristic equation exactly 
where specified. The steady-state error specification is satisfied 
by determining the relationship of the forward gain (Kt) to other 
system parameters. The bandwidth is found from the equation 

T 

Q ^ Kt TnHd(jCo) Cn(.o)) (3.22) 

R(jCo) 

Thus the bandwidth of the system is a function of the location of 
the roots, but the steady-state error solution is not. By judicious 
choice of his root locations, the designer may be able to meet both 
error and bandwidth specifications, or to determine whether a band- 
width specification conflicts with a stability or root dominance 
criterion. 
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The process of finding the proper root locations is Iterative, 
laborious, and time-consuming , requiring a complete solution of the 
system equations after each movement of the roots. For higher order 
systems, the assistance of a digital computer may be necessary. 
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4. Limitations of the Method 

Steady-State Error for a Ramp Input 
Steady-state error, as defined in Section 3, is the limit as 
s approaches zero of the product of s and the error function 
Error (s), where 



Error(s) 



v[ TdHd(6)tCd(s) + KtKfTnHn(s)*Cn(s) - KtTnHd(s)»Cn(sf | 
lnput^s;|^ TdHd(s)‘Cd(s) + KtKf TnHn(s)Cn(s) J 



For a unit ramp input, i.e. constant velocity signal of unity 

2 

magnitude. Input (s) « 1/s , Steady-state error becomes 




This limit approaches infinity as s approaches zero, unless an s can 
be factored from the numerator polynomial of the term in brackets, i.e. 



TdHd * Cd + TnHn • KtKfCn 
o o o o 



TnHd « KtCn 
o o 



(4.2) 



Equation 4.2 may be set equal to zero in several ways: 

1) The possibility of requiring that Cn^ or Cd^ be zero must 
be disregarded, as this would fix the value of a variable 
required for solution of the system of simultaneous equa- 
tions. 

2) Hn^ may be required to be zero, i.e. tachometer feedback, 
but the steady-state error expression resulting from this 
procedure introduces two new variables KtCn^ and KtCn^^ 

in a single equation 

Ess = f(Cd , Cd, , KtKxCn , KtCn , KtCn, ) 

' o’ 1’ o* o’ 1"^ 

and solution for Independent Kt is Impossible. 



Kt = 



Kt Cn 

c 

Cn 



Kt = 



Kt Cnj 
Cn, 
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3) Td^ may be required to be zero, i.e. a pole at the origin 

of the open loop transfer function. Then 

Tn Hn KtKfCn - Tn Hd KtCn = 0 

o o o o o o 

and _ Hd^ ; i.e. unity feedback. 

HiT 

o 

If this procedure is used Kf is no longer an independent 
variable and solution is impossible. 

4) If Hd^ is required to be zero, 

Tn Hn KtKfCn = 0 

o o o 

or Tn must also be zero, and an s factors from the denomln- 
o 

ator as well as from the numerator of Error (s). 

To date, the author has found no expression for steady-state 
error to a ramp input which is compatible with the system of 
simultaneous equations previously derived. 

Multi-path Feedback 




Figure 4 

Control System with Two Feedback Paths 
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The transfer function for the control system illustrated in 



Figure 4 may be written 



F(s) 



Kt T(s)«C(s) 

1 + Kt T(s)«C(s) [Kf H(s) + Kg] 



which reduces to 



(4.3) 



. . _ Kt Tn Hd(s)*Cn(s) .. .. 

" TdHd(s)»Cd(s) + KtKf TnHn(s).Cn(s) + KtKg TnHd(s)#Cn(s) 

The system of simultaneous equations obtained by equating coefficients 

of the characteristic equation will be linear in the variables (Cd^, 

Cdj^, — Cd^), (KtKfCn^, KtKfCn^, — KtKf) and (KtKgCn^, KtKgCn^^, 

KtKg) . The equations may be solved to obtain values for the 

variables composed of grouped parameters. However, for this system, 

and other multi-path feedback systems of various configurations, it 

has not proved possible to solve for independent values of the 

Independent parameters. . For example, in the system above KtKfCn^, 

KtKf, KtKgCn^, and KtKg are independent variables. But the value of 

Cn obtained from 
o 

Cn » KtKfCn 

o o 

KtKf 

is not necessarily the value of Cn^ obtained from 

Cn = KtKgCn 

o o 

KtKg 

Systems of this type have defied solution by this method. 
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5. Sample Solutions 
Example 5.1 



(y\- 


Kt 




1 








(s+l) ( 8 - 2 ) ( 8 + 8 ) (s+15) 




0(8) 




Figure 5.1 

Block Diagram for Example 5.1 

The open loop transfer function given above is unstable due 
to the pole in the right half plane (s » 2). Design a filter 
compensator and find values for forward and feedback gain which 
will make the system shown in Figure 5.1 absolutely stable. 

Dominant roots are specified by 5 “ 0.707 and oj^ ® 2.5. 

No additional steady-state error or bandwidth requirements 
are specified. 
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NOTE 



INPUl DATA 

«**<•****•»*** 



OPt-N LOrp TRANSFER FUNCTION - T(S) 

NUMERATOR OF OPEN LOOP 
TRANSFER FUNCTION - TIM(S) 

TN(S)= l.*'CCE + ^’0 

DENOMINATOR GE OPEN LOOP 
TRANSFER EUNCTION - TD(S) 

MJMGER OF ROCTS = A 

ROOT M(;. 12 3 4 

REAL PT. -.lOL+01 .2GE + :-l -.aOr + ':i -.15E+02 

IMAG PT. .ROF + :C .,'OE + :':' .COi. + CO ."CE+.H' 

COEFFICIENTS OF EXPANDED POLYNOMIAt 
TD(S)= l.:'C0F+:0 si 4) 

+ 2.20CCi+ 1 S( 3) 

+ 9.5CK0E + CI S( 2) 

+ -1.66CCi: + <‘2 SI 1) 

+~2.4 .’CwF + .2 

- PROGRAM HAS READ POLE(S) 

IN THE RIGHT HALF OF THE 
COMPLEX plane. 
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FEEDBACK PATH TRANSFER FUNCTION - H(S) 



X 



(0) ORDER FEEDBACK IS SPECIFIED. 

NUMERATOR OF FEEDBACK PATH 
TRANSFER FUNCTION - HN(S) 

HN(S)= l.OOOOE+00 

DENOMINATOR OF FEEDBACK PATH 
TRANSFER FUNCTION - HD(S) 

HD(S)= l.OOOOE+CO 



CLOSED LOOP TRANSFER FUNCTION 
G(S) - UNCOMPENSATED 

NUMERATOR CF CLOSED LOOP 
TRANSFER FUNCTION - GN(S) 

GN(S)= KT* l.GOOOE+00 

DENOMINATOR OF CLOSED LOOP 
TRANSFER FUNCTION - GDIS) 

GD(S)= l.OOOOF+rO S( A) 

+ 2.2000E+01 S( 3) 

+ 9.50C0E+01 S( 2) 

+-1.6600E+02 S( 1 ) 

+ (-2.A000E+02 + KT*KX» l.OOOOE+00) 
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STEADY STATE -kkU!< SPECIFICATION 



NO STEADY STAT.- ERkOR REQUIREMENT 
IS SPECIFIED. 



BANDWIDTH SPECIFICATION 

NO BANDWIDTH RFDUIREMENT IS SOECIFIED. 

NOTF - data INDICATES NO BANDWIDTH OR PARAMETER 

DEPENDENT STEADY STATE ERROR SPECIFICATIONS. 
UNDER THESE CONDITIONS THE PROGRAM WILL NOT 
SOLVE FOR FORWARD GAIN AND FEEDBACK GAIN 
SEPARATELY. SOLUTION IS TOR THE PRODUCT OF 
THE GAINS. (KT*KX) 



ROUTS SPECIFIED FOR 
CLOSED LOOP SOLUTION 



ZETA = .7G7CL + i'U OMEGA = 2.5GO E+00 

NUMBER OF ROOTS = 2 

ROOT NO. 1 2 

REAL PT. -.1BE + - 1 -.18E+''! 

IMAG PT. .lOh + Ol -.18E + ''-! 
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SOLUTION OF THE PROBLEM 



PROGRAM REQUIRES 7 ROOTS FOR SOLUTION. 



SOLUTION OF THE SYSTEM OF LINEAR EQUATIONS 



VARIABLE 






PARAMETERS 


RELATIVE ERROR 


VARIABLE 


( 


1 ) 


CO( 0) 


7.358E-09 


PERCENT 


VARIABLE 


( 


2) 


CD( 1) 


.OOOE-01 


PERCENT 


VARIABLE 


( 


3) 


CD( 2) 


.OOOE-01 


PERCENT 


VARIABLE 


( 


4) 


CD(3) 


-4.180E-09 


PERCENT 


VARIABLE 


( 


5) 


KT»KX*CN (0) 


-1.032E-07 


PERCENT 


VARIABLE 


( 


6) 


KT*KX*CN(1) 


3.466E-07 


PERCENT 


VARIABLE 


{ 


7) 


KT*KX*CN (21 


1.393E-06 


PERCENT 


VARIABLE 


i 


8 ) 


KT*KX*CN(3) 


2.421E-06 


PERCENT 


- RELAT 


IVE 


: ERROR 


STATED ABOVE 







MEASURES THE ACCURACY OF THE 
SOLUTION FOR EACH VARIABLE. 
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SOLUTION FOR INDIVIDUAL PARAMETERS 



NUMERATOR OF CASCADE COMPENSATION 
CN(S)= i.o:<ooF+c': s( 3i 

+ 1.1H67E+01 SI 2) 

+ 3.'«t32GF+?l SI 1) 

2.7105E + O1 

ZEROS OF THE FILTER 

ROOT NO. 1 2 3 

REAL PT. -.BOf+01 -.13E+01 -.25E+01 

IMAG PT. .COE+nC .OOF+3D .OOE+00 

NOTE - ACCURACY OF ROUT FINDING SUBROUTINE 
HAS BEEN CHECKED BY EXPANDING THESE 
ROOTS, AND COMPARING THE COEFFICIENTS 
OBTAINED WITH THE COEFFICIENTS OF THE 
ORIGINAL POLYNOMIAL. 

ERROR NOTED IS LESS THAN ONE PERCENT 
FOR ALL COEFFICIENTS. SOLUTION FOR 
ROUTS IS CONSIDERED ACCURATE. 
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DENOMINATOR OF CASCADE COMPENSATION 



CO(S)s l.OOOOE+00 S( 31 
+ 1.5117E+01 S( 2) 
1.4418F+02 SI 1) 
+ 2.9161E+02 



POLES OF THE FILTER 



ROOT NO. 


1 


2 


3 


REAL PT. 


-.63E+01 


-.63E+01 


-.26E+01 


IM^G l»T. 


-.85E+01 


.85E+01 


.OOE+00 











NOTE -? ACCURACY HAS BEEN CHECKED t AS BEFORE. 

ERROR NOTED IS LESS THAN ONE PERCENT 
FOk ALL COEFFICIENTS. SOLUTION FOR ' 
IS CONSIDERED ACCURATE. 

* - y 

PRODUCT OF GAINS IKT*KX) » 5.4508E^03 

/• . 

FORWARD' GAIN (KT) IS 5.4508E+03t 
WITH "UNITY FEEDBACK GAIN. (KX=1.0I 



I 
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CLOSED LOOP TRANSFER FUNCTION 
OF SOLVED PROBLEM 

NUMERAIOR POLYNOMIAL 
GN( S )= 5.A5G8F + 03 S( 3 ) 

+ 6.A686E+CA S( 2) 

+ 1.87C7t + C5 S( 1 ) 

+ 1.A775L+C5 

ZEROS OF THE TRANSFER FUNCTION 
ROOT NO. 1 2 3 

REAL PT. -.60E+01 -.13E+n -.25F+:i 

IMAG PT. .CCr+:0 .OCE+O: .OOL+00 

NOTE - ACCURACY FAS BEEN CHECKED, AS BEFORE. 

ERROR NOTED IS LESS THAN CNE PERCENT 
FOR ALL COEFFICIENTS. SOLUTION FOR 
ROOTS IS CONSIDERED ACCURATE. 
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DENOMINATOR POLYNOMIAL 



GD( S)= 


l.OOCOF+CO 


S( 


7) 








•f 


3.7118E+01 


S( 


6 ) 










5.7177F+C2 


S( 


3 ) 










4.7338E+C3 


S( 


A ) 






' 




2.281AE+0A 


S( 


3) 










6.A826F+C4 


S( 


2 ) 








+ 


1.3A06F+05 


S( 


1 ) 










7.7761E+C4 












POLES OF THE SYSTEM AND ROOTS 
CF THE CHARACTERISTIC EQUATION 








ROOT NO. 


1 




2 


3 


4 


5 


REAL PT. 


-. ISE + Ol 


- 


.18E+01 - 


. 53E+01 


-.53E+D1 


-.71E+01 


IMAG PT. 


. 18E+01 


- 


.18E+01 


. 18E+00 


-.18E+00 


.18E+00 


ROOT NO. 


6 




7 








REAL PT. 


-.71E+01 


- 


.88F+01 








IMAG PT. 


-. 18F + :0 




. 0 0 E + 0 3 
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Example 3.2 




Figure 5.2 

Block Diagram for Example 5.2 

Cascade compensation and the tachometer feedback device given 
are to be used in the closed loop system shown in Figure 5.2. The 
maximum permissible steadystate error for a step input is to be 
10 percent. 

Dominant roots are specified by ^ 0.5 and oj^ = 6. 

No bandwidth requirement is specified. 
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INPUT DATA 



OPEN LOOP TRANSFER FUNCTION - T(S) 

NUMERATOR OF OPEN LOOP 
TRANSFER FUNCTION - TN(S) 

TN(S)» l.OOOOE+00 S( 1) 

+ 5.0000E+00 

DENOMINATOR OF OPEN LOOP 
TRANSFER FUNCTION - TO(S) 

TO«SI» l.OOOOe+00 S< 4I 

♦ 3.3000E+01 S( 3) 

♦ 3.3600E+02 S( 2) 

+ 1.8800E+03 S( 1 ) 

+ 7.2000E+03 

RODfS OF TOTS) 

ROOT NO. 1 2 3 

REAt PT. -.90E+01 -.20E+02 -.20E+01 

IMAG PT. .OOE-i-00 .OOE-^OO -.60E+01 



A 

-.20E+01 

.60E<«-01 



FEFCPACK PATH TRANSFER FUNCTION - H(S) 



(1) ORDER FEEDeACK IS SPECIFIED. 

NUMERATOR OF FEEDBACK PATH 
TRANSFER FUNCTION - HN(S) 

FN(S)= l.OOCCfE + Oe S( 2) 

+ 7.GOCOi: + CC S( 1) 

+ 

DENOMINATOR OF FEEDBACK PATH 
TRANSFER FUNCTION - HDIS) 

HD(S)= l.C-!:'C0h: + ''0 S( 2) 

+ 3.*K'>COP4ri S( 1) 

+ 2.16005+02 
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CLOSEn LOOP TPANSFtK FUNCTION 
G(S) - UNCCMPENSATFD 

NUMERATOR OF CLOSEO LOOP 
TRANSFER FUNCTION - GN(S) 

GN(S)= KT#( 1.000'>f+00 S{ 3) 

+ 3.530 'E + 01 SI ? ) 

+ 3.660CE+02 SI 1) 

+ l.'"80C'E + 03 ) 

DENOMINATOR OF CLOSED LOUP 
TRANSFER FUNCTION - GDIS) 

Gois)= 1 . . :-''cE+o: si 6) 

+ 6.3:0CC+01 SI 5) 



5A20>^ + ':3 SI 


4 ) 








1.9C88E+34 + 


KT»KX* 


1 .OOjOE+' ') 


SI 


3) 


l,3618F+''5 + 


KT*KX* 


1. 2000F + - 1 ) 


SI 


2) 


6.22;PF+.5 + 


KT*KX« 


3 . 50-30E+ : i ) 


SI 


1) 


1 .5552E+ •‘6 + 


KT»KX* 


,oo;oF-:i ) 
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STIrADY STATE ERROR SPECIFICATION 



ERROR SPECIFIED (ERROR / INPUT) IS .1 



BANDWIDTH SPECIFICATION 
NO BANDWIDTH REQUIREMENT IS SPECIFIED. 



Ri‘( TS SPECIFIED FOR 
Cl.OStC LOOP SOLUTION 



ZETA = .5'0CE + (0 



NUMBER OF 


ROIIS = 


ROOT NO. 


1 


REAL PT. 


-. 3 jE+ ' 1 


IMAG PT. 


.^2E+‘ 1 



OMEGA = 6.:'00' 

2 

.30F+01 

,52E+:i 



'OOE + C 0 



C+'J ! 
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SOLUTION OF THE PROBLEM 



( 

v'> 



C 

( 



( 

i 

( 

( 

( 

( 

C - 

(, 

( 

( 

( 

( 

r 

( 

c 

( 



PROGRAM REQUIRES 11 ROOTS FOR SOLUTION. 



STEADY STATE ERROR SOLUTION 

STEADY STATE ERROR IS .lOOOE+00 

RELATIVE SPECIFICATION 
ERROR IS 3.727E-06 PERCENT 

ERROR SPECIFICATION IS SATISFIED. 
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SOLUTIGN fiF THt SYSTEM OF LINEAR COUATIONS 



NOTE 



VAR I ABLE 




PARAMETERS 


RELATIVE ERROR 


VAR I ABLE 


( 1) 


CD( C ) 


.nc.)E-oi 


PERCENT 


VAR I ABLE 


( .) 


COI 1 ) 


-1 .A60E-09 


PERCENT 


VAR lABLE 


( 3 ) 


CD( 2) 


-^.589E-09 


PERCENT 


VARIABLE 


( A) 


CD( 3) 


) .759E-09 


PERCENT 


VAR lABLE 


( ) 


CD( A) 


.c'r E-01 


PERCENT 


VAR I ABLE 


( h) 


CD( 5) 


.G0‘'E-01 


PERCENT 


VARIABLE 


( 7 ) 


KT*KX*CN (0) 


-S .A2 .’E-09 


PERCENT 


VAR I ABLE 


( 8 ) 


KT*KX*CN ( 1 ) 


-1 .27 E-08 


PERCENT 


VARIABLE 


( 9) 


KT»KX*CN (2) 


-3.431E-08 


PERCENT 


VAR lABLE 


( l.t) 


KT*KX*CN(3) 


9.962F-08 


PERCENT 


VAR lABLE 


(11) 


KT*KX*CN (4) 


f .72 )E-07 


PERCENT 


VAR lABLE 


( 12 ) 


KT*KX*CN(5) 


3.H3E-06 


PERCENT 


VARIABLE 


( 13) 


KT»CN( 0) 


^.727E-06 


PERCENT 


- RELATIVE ERROR STATED ABOVE 
MEASURES THE ACCURACY OF THE 
SOLUTION FOR EACH VARIABLE. 
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SOLUTION FOR INDIVIDUAL PARAMETERS 



NUMERATOR CF CASCADE COMPENSATION 
CN(S)= l.OOCOt+OC S( 5) 

+ 6.7668E+01 S( 4) 

+ 1.7775E+03 S( 3) 

+ 2.2639E+04 S( 2) 

+ 1.3997E+05 S( 1) 

+ 3.3682t+C5 

ZEROS OF THE FILTER 

ROOT NO. 1 2 3 4 5 

REAL PT. -.20L+:2 -.12E+02 -.18C+32 -.B4E+01 -.92E+01 

IMAG PT. .COE+OO .OOE+00 .OOE+00 .OOE+00 .OOE+00 



NOTE - ACCURACY OF ROOT FINDING SUBROUTINE 
HAS BEEN CHECKED BY EXPANDING THESE 
ROOTS, AND COMPARING THE COEFFICIENTS 
OBTAINED WITH THE COEFFICIENTS OF THE 
ORIGINAL POLYNOMIAL. 

ERROR NOTED IS LESS THAN ONE PERCENT 
FOR ALL COEFFICIENTS. SOLUTION FOR 
ROOTS IS CONSIDERED ACCURATE. 
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DENOMINATOR OF CASCADE COMPENSATION 



NOTE 



CD(S)= l.OOCOE + 00 S( 5 ) 

+ 7.2‘'0Df + 01 S( A) 

+ 2. 591r‘ + 03 S( 3 ) 

+ 2.7312t-+H SI 2) 

+ 1.712AE+05 SI 1) 

+ 4.U^3Ar + j5 

POLES OF THE FILTER 

ROOT NO. 1 2 3 A 5 

REAL PT. -.2AL+,'2 -.2AE+02 -.89E + C1 -.78E + 01 -.78C + C1 

IMAG PT. -.llE + 02 .llE + :'2 .GCE+.O -.26E+01 .26t + jl 

- ACCURACY HAS BEEN CHECKED, AS BEFORE. 

ERROR NOTED IS LESS THAN ONE PERCENT 
FOR ALL CCEFF ICIFNTS . SOLUTION FOR 
ROOTS IS considered ACCURATE. 

FORWARD GAIN (KT) - 1.6020E+03 
FEEDBACK GAIN IKX) = 1.2011F+CQ 
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CLOStD LOOP TRANSFER FUNCTION 
f F SOLVED PROBLEM 

NUMERATOR POLYNOMIAL 



GN(S)= 


1.A0PGE+C3 


SI 


3 ) 




1 .6AA7E + C5 


SI 


7) 




7.2280L+06 


SI 


6 ) 




1.7734F+r8 


SI 


5 ) 




2.6528F+C9 


SI 


A ) 




2.A737E+10 


SI 


3 ) 




1.AI/12F+11 


SI 


2) 


+ 


^.3966E+] 1 


SI 


I ) 




5.8275F+1 1 







ZEROS OF THE TRANSFER FUNCTION 



ROOT 


NO. 


1 


2 


3 


A 


5 


REAL 


PT. 


- . 20E + 02 


-.18E+02 


-. 18C+02 


-. 12E+02 


-.12E+02 


IMAG 


PT . 


.COE+^0 


-.46E-)2 


.A6E-02 


-. 36E-02 


.36E-02 


ROOT 


NO. 


ft 


7 


8 






REAL 


PT. 


-.8AE+01 


-.92E+01 


-.50E+01 






IMAG 


PT. 


.OOE+00 


.OOE+03 


. OOE + C0 







NOTE - ACCURACY HAS BEEN CHECKED, AS BEFORE. 

ERROR NOTED IS LESS THAN ONE PERCENT 
FOR ALL COEFFICIENTS. SOLUTION FOR 
ROOTS IS CONSIDERED ACCURATE. 
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DENOMINATOR I'ULYNOMAL 



GDI S)* 


i.'joc :e + ''0 


S( 11 ) 




1.35Cwr+0? 


S( 10 ) 




e.l371E+C3 


S( 


9 ) 


+ 


2.89C7E+05 


S( 


8 ) 




6. 73( dE + E 6 


S( 


7 ) 


+ 


1 . ; , H 1 0 r + ■' 8 


S( 


6) 


+ 


1.2275E+C9 


Si 


5 ) 




9.93A5E+C9 


Si 


A ) 




5.b663l +10 


Si 


3) 


-I- 


2.2290F+11 


Si 


2 ) 




5. A 7991 +1 1 


Si 


1 ) 




6.A 750F+ 1 1 







POLES OF The system anij roots 
CF THE CHARACTFRlSriC EQUATION 



ROOT 


NO. 


1 


2 


3 


A 


5 


REAL 


PT . 


-.30E+01 


- ,3^E+ j1 


-.90l + 'U 


-.9CE+01 


12E+02 


IMAG 


PT. 


.52E + '*'l 


-.52E + 71 


.52E+CO 


-.52E+00 


.52E+00 


ROOT 


NO. 


6 


7 


8 


9 


i:. 


REAL 


PT. 


-. 12E + E 


-.15E+J2 


-. 15E+02 


-. l'8E+02 


-.18E+C2 


IMAG 


PT. 


- .52L + J ' 


.52E+00 


-.52L+70 


.52E+00 


-.52E+00 


ROOT 


NO. 


11 










REAL 


PT. 


-.21E+02 










IMAG 


PT. 


.OoE+CO 











51 



Example 5.3 




Figure 5,3 

Block Diagram for Example 5,3 

The gain of the open loop device shown in Figure 5.3 is less 
than -20 db for all frequencies. Design a filter compensator and 
find values of forward and feedback gain which will provide positive 
gain at low frequencies and a bandwidth of zero db at 30 radians per 
second. 

Dominant roots are specified by ? 0,5 and “ 10 • 

No steady-state error requirement is specified. 
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[•\PUr DATA 

«»««««'«««»«« 



OPfcN U f TPANISFER FUNCTIUM - r(S) 

NUMFRATLiR C F LOOP 

TRANSFER FUrK TION - TN(S) 

TN(S)= 6.-> Cj.-+;1 



rENUMi';ATC.l (r i.PF) LOtJP 
TRANSFER FU FTION - TD(S) 



TD( s)= + v' S( ^ ) 

+ 3. ■ ■ 1 S{ n 

+ 2. n. •: + 2 S( 2 ) 

+ 7.<^- ^ JL + . ? S( 1 ) 

<■ S . / o ' + .:’2 

ROUTS t F TD( S) 

A 

25E + S1 
.OOE+OO 



ROUT 


NO. 


1 


7 


3 


WEAL 


PT . 


- . 12.- + S 2 


- . 1 5F + ':-2 


lAF + Sl 


IPAG 


PT . 


Ct +:■: 


. )0E + ::. 


.COF+00 
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FEEDBACK PATH TRANSFER FUNCTION - H(S) 



(0) ORDER FEEDBACK IS SPECIFIED. 

NUMERATOR OF FEEDBACK PATH 
TRANSFER FUNCTION - HN(S) 

HN(S)= l.OOCOE+00 

DENOMINATOR OF FEEDBACK PATH 
TRANSFER FUNCTION - HD(S) 

HD(S)= l.OOCOE+00 S( 1 ) 

+ 5.00C0E+C0 
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CLOSrO LUOP TRANSFER FUNCTION 



G(S) - UNCOMPENSATED 
NUMERATOR CF CLOSEu LOOP 



TRANSFER FUNCIION - 


■ GN(S) 


GN( S )= KT«( 6.3C0 


<F + 01 S( 1) 


+ 3.150- 


,'E+02 ) 



DENOMINATOR (F CLOSTD LOOP 



TRANSFER FUNCTION - 


- GDI S ) 


GDI S )= 1 J 


SI 5) 


+ 3 . 5 ■' V. F + L 1 


SI 4) 


+ A.23!‘:-+'? 


SI 3 ) 


+ 2.i:5Ci+ 3 


SI 2 ) 


+ A.?7f0f + :'3 


SI 1 ) 



+ ( 2.eS.rt-+>3 + KT*KX* 6.30:-OE + 'l) 
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STEADY STATE tKRG.{ SPEC I F IC AT I ') 



NO STEADY STATE ERROR REQUIREMENT 
IS SPECIFIED. 



BANDWIDTH SPECIFICATION 

GAIN OF .OOOE-Ol OB AT FREQUENCY 
OF 3.000E+01 RAD/SEC IS SPECIFIED. 



ROOTS SPECIFIED FOR 
CLOSED LOiiP SOLUTION 



ZETA = .5000E+00 OMEGA = 

NUMBER OF ROOTS = 2 
ROOT NO. 1 2 

REAL PT. -.50E+01 -.50E+01 

IMAG PT. .87E+01 -.87E+01 



i.ooo:e+oi 
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SULUllON OF THE PROBLEM 



PROGRAM RECUIRiS 9 ROOTS FOR SOLUTION. 

BANDWIDTH SOLUTION 

GAIN IS .iiCvI-Ol DB AT FREQUENCY 3.‘OOE + Ol RAD/SEC 
RELATIVE ERROR IS .OOCE-JI PERCErJT 
BANm. IDTF SPFC If ICM ION IS SATISFIED. 
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NOTE 



1 FML SYSTEM OF LINbAR EQUATIONS 
PARAMETERS RELATIVE ERROR 



VAr. ' 


( 


1 ) 


CO( 0) 


-2.355E-09 


PERCENT 


VAR I ABLE 


( 


2 ) 


CD( 1 ) 


-5.282F-09 


PERCENT 


VAR I ABLF 


( 


3) 


CD! 2) 


-7.561E-08 


PERCENT 


VAR I ABLfc 


( 


A ) 


CU( 3) 


-3.551E-07 


PERCENT 


VARIABLE 


( 


S ) 


CD( A) 


1.5S9E-06 


PERCENT 


VAR I ABLE 


( 


t) ) 


KT*KX*CN(0) 


7.65AE-05 


PERCENT 


VAR I ABLE 


( 


7) 


KT*KX*CN( 1) 


1 68E-03 


PERCENT 


VARIABLE 


( 


8 ) 


KT*KX*CN( 2) 


8.962E-03 


PERCENT 


VAR lABLF 


( 


■7) 


KT*KX*CN(3) 


S.280E-02 


PERCENT 


VARIABLE 


( 1 


J) 


KT*KX*CN ( A) 


2.AA1E-01 


PERCENT 


- RELAT 


IVE 


ERROR 


STATED ABOVE 







MEASURES THE ACCURACY OF THE 
SOLUTION FOR EACH VARIABLE. 
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SI LIjTinN FOK INDIVIDUAL PARAMLTEKS 



MJMFRATUR i,F CASCADE COMPENSATION 
CN( S )= 1. iT'-i + j. s( A ) 

+ 3. . 7’-+ 'l S( 3) 

+ 5.A/Mr+.? S( 2) 

+ 3 .A 0 +'• 3 S( 1 ) 

+ « . A 9 5 + : 3 

ZEROS IN- THE Fill hC 

ROOT NO. 1234 

»FAL PT. -.12<+'2 -.1dE+ 2 -.6LE+C1 -.61E+ 1 

IMAG PT. + 3 .03E+'"' -.36L + 31 .36E+ul 

NOTl - ACCUI^ACY (» ■•(•'I FliJDIN', SUBROUTINE 

HAS HEE\ CHECKED BY EXPANDING THESE 
ROOTS, AND cnrPAKiNG THr COEF F I C I FIN T S 
OBTAINED WITH THE COEFFICIFNTS OF THE 
ORIGINAL POLYNOMIAL. 

FPRi K NGThO IS LESS THAN ONE PERCENT • 

FOR ALL CCLi FICIENTS. SULUTIOrj FOR 
ROOTS IS rO.SIDEREO ACCURATE. 
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DENOMINATOR OF CASCADE COMPENSATION 



NOTE 



CO(S»= 1.002AL+00 S( A) 

+ 1.25C0E<-02 S( 3) 

+ 6.35A3F>03 S( 2» 

+ 1.6789E+05 S( 1 ) 

+ 2.AA27E+06 

POLES OF THE FILTER 

ROOT NO. 1 2 3 

REAL PT. -.A7F.+02 -.A7E+j 2 -.15h+32 

IMAG PT. -.17F + 02 .17E+02 -.27^^ + 02 

- ACCURACY HAS BEEN CHECKED, AS BEFORE. 
ERROR NCTEO IS LESS THAN ONE PERCENT 
FOR ALL COEFFICIENTS. SOLUTION FOR 
ROOTS IS CONSIDERED ACCURATE. 

FORWARD GAIN (KT) = A.2A11E+04 
FEEDBACK GAIN (KX) = 7.1267E+00 



A 

. 15E+02 
.27E+D2 
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CLDSEl) LOliP TRA.NSFER FUNCTION 
I F SULVLD PROBLEM 



NUMERATOR POLYfiOMIAL 
CN( S )= 2.6719!. +‘ 6 S( 5 ) 

+ 1 . li^25!'+08 S( AJ 
+ 1.96-^ r + .. 9 S( 3) 

+ 1 .633or + lC- S( 2 ) 

+ < . 8 1 1 + 1 : s ( 1 ) 

+ i.i4;'2--m 

ZEROS (F FH= TRANSFER FUNCTION 



ROOT 


NU . 


1 


c. 


3 


4 


!:> 


REAL 


PT . 


- . 6 1 r + 0 1 


-.6 lE + Dl 


-. 15E+G2 


-. 12E + 02 


-•50rr^Gl 


I NAG 


PT. 


- .36r +' 1 


.36E+01 


.0C(-+ :o 


.0OE+C‘0 


.0 JL+00 



NOTE - ACCURACv HAS BEEN CHECKEOtAS HEFORF, 

ERROR NOTED IS GREATER THAN ONE PERCENT 
FOR AT LEAST ONL COEFFICIENT. SOLUTION 
FUR ROOTS IS NOT CONSIDERED ACCURATE. 
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DENOMINATOR POLYNOMIAL 



GDI S) 


ZZ 


1.C0C0L- + 00 


SI 


9) 












1.60C0E+02 


SI 


8) 












1.1152t+0A 


SI 


7 ) 












A.A527E+C5 


SI 


6 ) 












1. 1274E + 07 


SI 


5 ) 










-•> 


1 .89A7E + 08 


SI 


4) 












2. 1A74E + 09 


SI 


3 ) 












1.62C1F+10 


SI 


2 ) 












7.5736F+10 


SI 


1 ) 












1.6983E+11 












POLES OF THE SYSTEM AND ROOTS 
CF THE CHARACTERISTIC EQUATION 








ROOT 


NO. 


1 




2 


3 


4 


5 


REAL 


PT. 


-.50E+01 


- 


.50E+01 - 


. 15E+02 


-. 15E+02 


-.20E+02 


IMAG 


PT. 


.e7E+Cl 


• 


.87E+01 


.87E+00 


-.87E+00 


.87E+C0 


ROOT 


NO. 


6 




7 


8 


9 




REAL 


PT. 


-.20E+02 


• 


.25E+02 - 


.25E+02 


-.30E+02 




IMAG 


PT. 


-.87E+0C 




.87E+00 - 


.87E+00 


.OOE+00 
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(t 



Example 5.4 
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00 
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o 








•X 


o 








o 
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cd 


CO 
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0) 


O 








<y 


P4 
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d 








d 


Cd 






<r 


cd 
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'O 






in 
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o 
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}4 
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d 
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(0 


cn 


d« 


cd 


u 


0> 
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o 


r- 
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^4 


44 


r4 


o 
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cO 


O 
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03 


• 


(0 
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o 
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44 






P 


C 


T) 


44 


u 
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f4 


d 
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o 


T) 


d 


•rl 
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43 


p 
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CJ 
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U 
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0) 


d 


<0 


CO 




44 
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<u 
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•rl 
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44 


te 


(d 
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CO 


1 
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>> 
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00 


a 
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•rl 


0) 


Cd 
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CO 


44 
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O 


CO 
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>> 
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CO 
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d) 
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44 
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INPUT DATA 



OPEN LOOP TRANSFER FUNCTION - T(S» 

NUMERATOR OF OPEN LOOP 
TRANSFER FUNCTION - TN(S) 

NUMBER OF POINTS = 2 

ROOT NO. 1 2 

REAL PT. -.6'L + Ol -,13E + ''2 

IMAG PT. .OSL + OO .OLE + .'J 

COEFFICIENTS OF EXPANDED PilLYNOMIAL 
TN(S)= 1.j:LuL+CO S( 2) 

+ 1.9jCOb+0l S( 1) 

+ T.eOCOE+Ol 
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CENUMINATOR (’F UPFN LOOP 
TRANSFER FU.CTICN - TU(S) 

NUMBER OF ROOTS = t, 



ROOT 


Nl<, 


1 




2 


3 


A 


5 


REAL 


PT . 


. V t. + C •; 


- 


.5C£+:i 


-.5;.E + 31 


-.9CE+01 


-.15L+V2 


I MAG 


PT . 






.30E+ J1 


30! +:i 


. 'OE+00 


.80E+C1 


ROOT 


NP. 














REAL 


PT . 














I MAG 


PT . 














COEFFICIENTS Of fX 


-'ANnEU POLYNOMIAL 






TD( S ) 


= 


1. r . Vh + sc 


S( 


6 ) 










+ 


^.gcco. + " 1 


S( 


5 ) 












9.8TC01 + OT 


S( 


A ) 










+ 


9.S17JL+03 


S( 


3 ) 










+ 


^ . 5 : 16L+ V.A 


S{ 


2 ) 










+ 


8. RA39L+0A 


S( 


1 ) 












.■’XCL-.)! 
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FEEDBACK PATH TRANSFER FUNCTION - H(S) 
(0) ORDER FEEDBACK IS SPECIFIED. 

NUMERATOR OF FEEDBACK PATH 
TRANSFER FUNCTION - HN(S) 

HN(S)= l.OOCOE+00 

DENOMINATOR OF FEEDBACK PATH 
TRANSFER FUNCTION - HD(S) 

HD(S)= l.OCOOE+00 
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CLOSEH LOl;P TRANSFER FUNCTION 
G(S) - UNCCMPENSATEU 

NUMERATOR OF CLOSED LOOP 
TRANSFER FUNCTION - GN(S) 

GN(S)= KT*( l,0r0.)E + 00 S( 2) 

l.PC: -f+Ol S( 1) 

+ 7.6SC. E + Ol ) 



CENOKIf.ATOR CF CLOSLD LOOP 
TRANSFER FUNCTION - GDIS) 

GD(S)= 1.:' :oF+'’.o s( 6) 

+ ^.9 . + .'1 SI 5) 

+ :t +:'2 SI 4) 

+ 9.9 170) +c. 3 SI 3 ) 

+ I A.5;16L + j4 + KT*KX» 1 . CC CE + ''.) SI 2) 
+ I 8.e4 24E+.A + KT*KX* 1.90.-0E+'l) ST 1) 



+ I .rO;.OL-31 + KT*KX« 7.80,v;E:+ 1) 
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STEADY STATE ERROR SPECIFICATION 



ERROR SPECIFIED (ERROR / INPUTI IS 



BANDWIDTH SPECIFICATION 

GAIN OF -3.0f0E+00 DB AT FREQUENCY 
OF 2.C-OOE+(''1 RAO/SEC IS SPECIFIED. 



ROOTS SPECIFIED FOR 
CLOSED LOUP SOLUTION 



ZETA 


~ • 


707LE+00 


OMEGA 


NUMBER OF 


ROOTS = 2 




ROOT 


NO. 


1 


2 


REAL 


PT. 


-.28E+01 


-.28E+01 


IMAG 


PT. 


.28L+01 


-.286+01 



.lOOOE+00 



OOOOE+00 



SOLUTICN UF THE PROBLEM 

• •»««*«*•»««»•» » »•« «»»«»• 



PROGRAM REGUIRfS 11 ROOTS FOR SOLUTION. 



BANDWIDTH SOLUTION 

ITERATIVE PRCCirDURE FOR BANDWIDTH SOLUTION 
PROCESS IS CONVERGING 

ITERATION GAIN-DB RELATIVE ERROR 

31 -2.64f+C0 1.20E+C1 PERCENT 

36 -3.--'2L+00 -6.65E-01 PERCENT 

GAIN IS -3.S2 jb + '-'C' OB AT FREQUENCY 2.'00E + 01 RAD/SEC 
RELATIVE ERROR IS -6.655E-01 PERCENT 
BANDWIDTH St^EC I F I C A TI ON IS SATISFIED. 

STEADY STATE ERROR SOLUTION 

STEADY STATE ERROR IS .lOOJE + 'JD 

RELATIVE SPECIFICATION 
ERROR IS 3.27AE-0R PERCENT 

ERROR SPECIFICATION IS SATISFIED. 



* 
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SOLUTION OF THE SYSTEM OF LINEAR EQUATIONS 



VARIABLE 

VARIABLE 


( 1) 


PARAMETERS 
CD( D) 


RELATIVE ERROR 

-2.290E-09 PERCENT 


VARIABLE 


( 2) 


C0( 1 ) 


.OOOE-01 


PERCENT 


VARIABLE 


( 3 ) 


C0( 2) 


.QOOE-01 


PERCENT 


VARIABLE 


( 4) 


CD( 3) 


1.689E-09 


PERCENT 


VARIABLE 


{ 5) 


CD( A) 


.OOOE-01 


PERCENT 


VAR I ABLE 


( 6J 


CD( 5) 


3.595E-09 


PERCENT 


VARIABLE 


( 7) 


KT*KX*CN(0) 


1 .032E-08 


PERCENT 


VARIABLE 


( 8) 


KT*KX*CN( 1) 


1 .033E-07 


PERCENT 


VARIABLE 


{ 9) 


KT*KX*CN (2) 


1.555E-06 


PERCENT 


VARIABLE 


( 10) 


KT*KX*CN( 3) 


1 .347E-05 


PERCENT 


VARIABLE 


( 11) 


KT*KX*CN(4) 


8.095E-05 


PERCENT 


VARIABLE 


{ 12) 


KT*KX*CN(5) 


4.292E-04 


PERCENT 


VAR I ABLE 


(13) 


KT*CN( 0) 


-5. 123E+00 


PERCENT 



NOTE - RELATIVE ERROR STATED ABOVE 
MEASURES THE ACCURACY OF THE 
SOLUTION FOR EACH VARIABLE. 
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SOLUTION FOR INDIVIDUAL PARAMETERS 



NOTE 



NUMERATOR OF CASCADE COMPENSATION 
CN(S)= l.O^COt + rJ S( 5) 

+ A.5730h+Ll S( 41 
+ 6.342AI.+ 02 S( 3) 

+ 6.9977r + '3 S( 2 ) 

+ 2.7‘>] 5t + .)4 S( 1 ) 

+ 4 • 7 7 9 8 f + 4 



Z^^POS OF THE FILTE.. 

ROOT NO. 1 2 3 4 

REAL PT. -.15F+:2 -.15E+^2 -.90E+j1 -.34E+j1 
IMAG PT. -.76E + ' 1 .76E+01 .OOL+CO -.27E+01 

- ACCURACY CF ROOT FINDING SUBROUTINE 
HAS BEEN CHECKED BY EXPANDING THESE 
ROOTS, AND COMPARING THE COEFFICIENTS 
OBTAINED MTH THt COEFFICIENTS OF THE 
ORIGINAL POLYNOMIAL. 

ERROR rjOTEl) IS LESS THAN ONE PERCENT 
FOR ALL CCEFF ICIENTS . SOLUTION FOR 
ROOTS IS COiSIDERED ACCURATE. 



5 

-.34E+01 
.27L- + 01 
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CENOMINATGR CF CASCADE COMPENSATION 



CD(S)= l.OoOOE+aO S( 5) 
+ 1.7057E+02 S( 41 
+ 1.1647t+04 S( 3) 
+ 4.C007E+05 S( 2) 
+ A.0045E+06 S( II 
+ 5.8576E+06 



POLES OF THF FILTER 
ROOT NO. 1 2 

REAL PT. -.80E+02 -.38E+02 

IMAG PT. .OCE+00 -.43E+02 



3 

-.38E+02 

.43E+02 



4 

-. 18E+01 
.OOE+00 



5 

-.13E+02 

.006+00 



NOTE - ACCURACY HAS BEEN CHECKEOtAS BEFORE. 

ERROR NOTED IS LESS THAN ONE PERCENT 
FOR ALL COEFFICIENTS. SOLUTION FOR 
ROOTS IS CONSIDERED ACCURATE. 



FORWARD GAIN (KT) = 2.6980E+C6 
FEEDBACK GAIN (KX) = l.llllE+00 
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CLOSl'D LOOP TKANSFER FUNCTION 
Tf SOLVED PROHLEM 

NUMERATOR PDLYDUMIAL 



GN( S ) 




2.6980L + 


S( 


7 ) 










+ 


1 .7464! +08 


S( 


6 ) 










+ 


4.6- SSC+Og 


S( 


5 ) 










4 


7.12691+10 


S( 


4 ) 










4 


6. ■ 8S3r + 1 1 


S( 


3 ) 










4 


3 . n 2 1 1. + 1 2 


S( 


2 ) 










4 


8.241 7! + 12 


S( 


1 ) 










4 


1 .0’: 59L + 1 3 












ZEROS 


OF 


THE TRANSFER 


FUNCTI CN 








ROOT 


NO. 


1 




2 


3 


4 


6 


REAL 


PT • 


-. 13l + 0? 


- 


.15E+02 - 


. 15E+02 


-.60E+01 


-.90t+01 


IMAG 


PT* 


. C 0 L + C 0 


- 


.76E+01 


.760+01 


.OOE+OC 


.OOE+00 


ROOT 


»\0 . 


6 




7 








REAL 


PI. 


“ .34c + ^^l 


- 


.34E+D1 








IMAG 


PT. 


-.2 71 + 1 1 




.27E+ )1 









NOTE - ACCURACY HAS HEEN CHECKED, AS EEFORE. 

ERROR NCTEO IS LESS THAN CNE PERCENT 
FOR ALL COEFFICIENTS. SOLUTION FOR 
ROOTS IS CU.NSID6REO ACCURATE. 
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DENOMINATOR POLYNOMIAL 



GDI S ) 




l.OOCOE+00 


SI 11 ) 










•I- 


2.1957E+02 


SIIO) 












2.0988F+0A 


SI 


9) 












1.1A80E+06 


SI 


8 ) 












3.9723t- + 07 


SI 


7) 












9.08C0F+08 


SI 


6) 












1.3910F+10 


SI 


5) 










+ 


1.A210E+11 


SI 


A ) 












9.A753t,+ ll 


SI 


3) 












3.96A6E+12 


SI 


2 ) 












9.67A3E+12 


SI 


1 ) 










+ 


1.1177E+ 13 












POLES OP THE SYSTEM AND ROOTS 
CF THE CHARACTERISTIC EQUATION 








ROOT 


NO. 


1 




2 


3 


A 


5 


REAL 


PT. 


-.28E+01 


- 


.28E+01 - 


.85E+01 


-.85E+01 


-. 17C+C2 


IMAG 


PT. 


.28E+01 


- 


.28E+01 


.28E+00 


-.28E+00 


.28E+00 


ROOT 


NO. 


6 




7 


8 


9 


10 


REAL 


PT. 


-.17E+02 


- 


.26E+02 - 


.26E+02 


-.3AE+02 


-.3AE+02 


IMAG 


PT. 


-.28E+00 




.28E+00 - 


.28E+00 


.28E+00 


-.28E+00 


ROOT 


NO. 


11 












REAL 


PT. 


-.A3E+02 












IMAG 


PT . 


.OOE+00 
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6. Conclusion 



The method presented has been shown to be fast and accurate, 
and capable of designing systems which will meet specifications 
for root dominance, steady-state error, and bandwidth. 

The principal disadvantage of the method presented, exclusive 
of the restrlctloEB mentioned In Section 4, Is the large number of 
filter sections required, particularly for higher order open loop 
systems. The filter sections tend to further complicate the 
originally complex systems. 

The possibility of reducing the number of filter sections 
employed In the final design Is worthy of further study. Approxi- 
mation of the characteristics of the cascade compensation may be 
accomplished with a filter of lower order, but analysis of the 
effects of such approximation on system specifications Is beyond 
the scope of this paper. 
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program COMPTECH is written in the fortran 60 LANGUAGE FOR USE 
WITH A CDC 1604 COMPUTER SYSTEM. THE TAPE LIBRARY AVAILABLE AT 
NAVAL POSTGRADUATE SCHOOL COMPUTER FACILITY PERMITS USE OF THE 
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FIRST FIELD - THE INTEGER ZERO (0). INDICATES THAT ROOTS 



SECOND FIELD - THE NUMBER OF ROOTS TO BE REW> FOR TN(S) 
MAXIMUM IS 9. 
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FIRST DATA CARO 

FOUR FIELDS OF 15 FORMAT 
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TEN SPACES* AND THE REMAINDER OF THE CARD IS LEFT BLANIC 
EACH CARD HAS SPACE FOR EIGHT COEFFICIENTS. ADDITIONAL 
COEFFICIENTS REQUIRE ADDITIONAL DATA CARDS. 



SAME FORMAT AS OAT A CAROS FOR TN<S)# WITH THE EXCEPTION 
THAT TOCS) MUST BE A RATIONALIZED POLYNOMIAL * I ,E* THE LAST 
COEFFICIENT READ MUST BE ONE (1*0>. 
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BANDWIDTH IS SPECIFIED AS A FREQUENCY IN RAOIANS PER SECOND AT 
WHICH THE SYSTEM WILL HAVE A DESIRED GAIN IN DECIBELS* 



EIGHTH DATA CARD 
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FORMAT (19X*13H ************ ) 

This section reads the numerator and denominator of the open loop 
transfer function, either as the roots or the coefficients of 
polynomials. 

PRINT 5 
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IF (MARK) 100 ♦100.9999 
100 CALL RTEXPD (NDTI .TDR.TDI *TD) 
IF (T0(0)) 130.110»130 
110 KTDO-0 

IF (TD(D) 130.120.130 



120 KTD1»0 

130 CALL DATAPR2 ( NOT I *TD*0*2 ) 

CALL RHPCHIC (0. 1*N0T I »TDR*KRHP#KOUT ) 
IF (KOUT-KRHP) 200 t 200 #9999 
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223 FORMAT (1H0*13X»37H THAN ORDER OF PRODUCT OF NUMERATORS, 
PRINT 2231, NNTI,LFPNH,NTPH 

2231 FORMAT ( IHO ,/ , 14X, 31H ORDER OF NUMERATOR PRODUCT IS ,I4,3H 



PRINT 2232* NDT I *NDHM1 *MTPH 

2232 FORNAT ( IHO*/ . 14X*33H ORDER OF DENOMINATOR PRODUCT IS *I2.3H 
1 I2*3H «= *12*// ) 

MARK»1 

224 IF (10-MTPH) 225*228*228 
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445 PRINT 450 

450 FORMAT ( 1H0»/// *10X» 28H DENOMINATOR OF CLOSED LOOP 
PRINT 451 

451 FORMAT (1H0.9X.27H TRANSFER FUNCTION - GDIS) ) 
CALL PLYMPY ( NDT I » TD*NDHM1 *HD.N1D»TDHD) 



6ALL PLYMPY { NNT I »TN » LFPNH »HN,N20* TNHN ) 
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NCT»0 

MCML=0 

AP»-1,0 

BLFJ=1.0 

CALL BANDSLV (N1N.TNHD»N1O,TDhD.N2D,TNHN,ECD0,EKxN0»EKN0» JRSPC 



1 FCK*FCJ»BASE»FAQ*RTRL*RTIM*B0MEG»GM6N*GMG0»GMGT»MARK) 
IF (MARfc) 773*773.9999 
773 IF <BDZ-GMGT) 785.766.775 
775 IF (GMGT-BOA) 820.766.766 
780 LTRY»11 
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911 PRINT 912 

912 FORMAT! 1H0*/*10X*42H BANDWIDTH SPECIFICATION IS NOT SATISFIED 
917 CALL ROOTGEN ( JRSPC* IMX*FCK.» FC J*8ASE *FAQ *RTRL» RT I M ) 

920 CALL RTEXPD ( I MX*RTRL *RT IM,BAR ) 

CALL 8ARARRG (I MX* BAR) 



CALL COPLUG (N10*T0HD»N2D.TNHN*NMAT .COLE) 

bXl«NMAT+l 

LX2-N1D+1 

€OLE(LX1»1)=ECDO 

60LE(LX1.LX2)»EKXN0 
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PRINT 9690 

9690 FORMAT (1H0,12X*28H SOLUTION FOR EACH VARIABLE. ) 

This section computes and prints out the gains and the 

COEFFICIENTS AND ROOTS OF THE FILTER. 
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SUBROUTINES DATAPR(0 TO 5) ARE PRINTOUT FORMATS FOR ROOTS AND 
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5 PRINT 10»NAB 

10 FORMAT (1H0*/*10X» 18H NUMBER OF ROOTS 
15 BO 20 I«1»NAB 



< 

I—* 

X 

i 



GO 

< 

z 

X 



CO 

fVJ R 

<M — 



O w 

CM U 

Z 



o 



CO I o 

< CO O CO 
Z < rH 

M Z M l~ 
CO -- CO Z 

< < 



fS:fS 







-1 








o 








a -- 








X 


X 






O CM 


CM 






UJ 


*> 






O CM 


CM 






z ^ 


♦— 1 






< 


•k 






a -- 








X </) 


1/} 






Ul 








X 


X 






U. CO 


CO 


O CM CM 




O 




i-H • • 


>• 


4 


4* 


♦-Hi 


X 


cO • 


• 


0\ r-H t— H 


i-H 


h- CM 


CM 


» LU UJ 




Z rH 


f-H 


00 O O 


I— 


liJ UJ 


UJ 


♦-H I-HI f— t 


0^ 


a 


Ql 


» » •>. 


o 


yj ^ 


fH 




u 


t-H 


•h 




") 


u. 




OCOh-COh- 




U. M 


II 


z < a < a 


rM 


UJ -- 




z z 


< 


O u> 


1/) 


>- r -1 X o 








o • < • < 


GO 


z 


a 


o <-H UJ I-H X A 


< 


X h- 


I— 


CV l-H t—4 HH yf\ 


Z 


>0 ctj 


CO 


R H 




CO < X 


< X 


X ^ X -• X o 




• o z r- 


z 



Qn #k 0*' 


Os >o 






CM 




X 


4 












•k — ^ 


•k 






Z A 




o 


* 




X 






X 


X X «-i 


X o 






0. -- -1 




^ o 


<n 


o 




GO 


o 


rg w CM ^ 


CM vO 






< A • 






A 


< 


I-H 




< 


I-H 


•♦Of »-H 


* 






CM A 






•k 


z 


* 




z 


* 


V fsi rg 








< w • 




V. 


o 


w 


X. 








» < • < 


• CO 






O rg A 






CM 


rsl 


* 




rsl 


* 




o < 






< 




o 




< 


o 




< 


o 


X X 


X z 






UJ ^ 




X 






X 






X 


rH • »— 1 • 








z z ►- 




rH 


CM 


•k 


iH 


A 


•k 


rH A 


w o -- o 


-- T 


CO 




o O O' 


o 


w 


1 


A 


w 


rH 


A 


w rH 


4 A 


CO < 




CM H- -• O 




>• CM 




I-H CO 


rH 


h— 1— 


1 - < z 




Z 3 U) U 




I— 


X 




h- 






h- 


^ 1— 


< z 


N 


rH 


O « O Z “» ►- < 


t-H 




< o »- 


< o 


X Z X z X - 


GO 


•-3 Q£ UJ — 


Z X 


w 


Z 


Z I- z 


X 


Of *-• (T 


0^ 


< 


R 


H- O 03 Z 


I-H 






I-H 


z 




I-H 


z 


O cr O Q£ o u. 


z 


-^OUJZD-'U.QCO 


u z 


o 


O Z 


o o 


Ii. a u. a 


u. ^ 


X 




A Z liJ U> CO CX. 


u. 


I-H 


0. 


u. 


CO Z 


U. lO 



X 

CM 

» 

CM 



(/) 



X 

CO 

CM 



UJ 

a. 



ii 



1/) 



z 

X 

GO 

< X 
m z r- 

• ^ X 

m 03 o 

lA ^ < 

• Z 

m w >v 
4^ Nl * 

< O 

-- X 

4- ♦ lA 
I O -- rH 
>• m i-H 



X 



-- z 

U. OC 

l-H OL 



< O 
X »- 
OC 

o o 

Ii. 



CO On O O O lA 

CM CM CO 4- A A 



O 

>o 



A O A O A 
«-l i-H CM CM 



O A O A O 

CO CO ^ 4" A 



103 



55 PRINT 60* AZ(NAB)*NAB 

60 FORMAT ( IHO */ . lOX* 7H HD(S)* *1PE12.A*3H S(*I2*2H) 
60 TO 115 



























CM 




























«k 




























CM 
























































«k 




X 




X 


X 


X 




















CM 




CM 


CM 


CM 
















l/> 




* 








» 




















CM 




CM 


CM 


CM 
















X 








►-< 


►— « 


















<n 




* 
























•k 












w 






X 








4- 


4- 




i/} 




1/) 


i/> 


cO 






CM 








• 


• 


















•k 








iH 


04 




X 




X 


X 


X 






CM 








rH 


rH 




CO 




cn 


cn 


cn 






♦-H 








LU 


lU 




* 




* 


* 








#k 








a. 


Q. 








4- 


4- 


4- 














rH 


rH 




• 




• 


• 


• 






CO 








•k 






CM 




CM 


CM 


CM 








•— * 












rH 




rH 


•H 


rH 






X 














UJ 




LU 


LU 


lU 






CO 




LU 










Q. 




a 


Q. 


Q. 






•k 




e> 




»- 


H- 




rH 




1—1 


rH 


iH 






4- 


4* 


X 




X 


X 




* 




m 


•k 


•k 






• 


• 


•k 
























rH 


rH 


LU 




it 


if 




u 




it 


H 


li 






rH 


rH 


-1 










«— 






— «► 


— 






LU 


LU 


CO 




lO 


LO 




i/> 




1/) 


lO 


lO 






a. 


Q. 


> 










w 






w 


W 






•H 


iH 


«k 




X 


X 




X 




o 


X 


o 






«k 


•k 


X 




o 


o 




u 




u 


00 e> 


00 O 


m 




-I- 


+ 


n 












GO 




m < 


< 


CM 








H-l O 




X 


X 




X 


< 


X 


O X X 


X X 


rH 




X 


X 


^ vO 




CM 


CM 




r- 


X 


c- 


rH • C- 


♦ c- 


•> 




CM 


CM 


•k 




fH 


rH 




* 




* 


m ^ » 


r-N. » 


o 






» 


CO O 




•k 


X • 




>< 




X 


m CO X 


(O X 


in 




• X 


X 


QJ: rH 


O 


X 


-) X 


o 


o 


CO 


o 


o < o 


< o 


iH 


rH 




^ vO 


Q. • 


4- 


o 


•H O 


<> 


rH 


< 


rH 


rH Z rH 


X -* 


•k 


X 


“> rH 


o 


< o 


«k 


rH 


•k rH 


* 




X 




♦ w ♦ 




o 


CD 


<w » 


w » 


h— rH 


O 


» 


lU » 


o 


N. 




N. 


m fsi N. 


rsl N. 


lA 


< 


tsl N. 


N4 N. 


< 


CM lU 


N. 


-1 


r- 


* 


N 


* 


o> < • 


< ♦ 


•H 


X 


< • 


< ^ 


D — ‘ 


• -J 


m 


GO • 


«k 


o 


< 


o 


o 


o 


rH 


•k 


o 


O 


X 


O CO 


O 


> o 


o 


X 




X 


— • X 


♦ X 


-- 1 


•H 


♦ X 


• X 


LU *0 


nj > 


X 


X 


r- 


rH m 


* 


rH in 


GO O rH in 


O *H 


rH CO 


N 


O rH 


O rH 


X -- 


•k 


rH O 


• rH o 





r-H iTl I O'-' 

rH 00 rH >- rH 

I- H- X H- 

XH-ZXI — 

cr ^ cc ^ (X 
oo<rooii-(ro 



O ► 

K r X 

-i (X, 

o cr o 

a. u. 



I < — 
00 X 
< » o 

Z rH ^ 

w X .-H 

GO 

U- < o 
X (Ji 



CO w liJ ^ 

— • rH O rH 

I H- X I- X 

CO Hr < Q£ 



o 

oc 

X^XX^Ct^C^CD 
M QGOOQGOUJXD 
")Q.U-U>aU.OGlUtO 



< X X X X ^ 

«-H 2 ^ 



I ^ o 

UJ u ^ 

e> n 

X ^ h- 



U. U. QG 

•ii HH a. 



^ rH O rH — 

rH in rH ^ 

h- h- 

o — 



< o »- < 
X H- X X 
QG — QC 

o o q: o 

U. ^ 0 . U- 



in 

r- 



o m 

00 00 



o in o 
o Ch o 



m o in in 
o iH rH nj 



o o o o 

CO 4* in 'O 



o o o 

rH CM CO 



o o 
^ m 



O U- 

ij) ^ 

o 

sO 



lOA 
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SO TO 110 
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CALL TOTSUM ( 12 * I »FACRL »FACI M, SPFAC ) 
J*12-I 



50 SP1(J)«SPFAC 
SP1(12)=1.0 
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10 L9«l 

CALL SUMPR9 ( L9 , I » NC »MAX » FACRL . FAC IM»SUMRL»SUMIM ) 
60 TO 110 



CALL SUMPR8 ( L8 » I * NC »MAX * FACRL ♦ FAC I M »SUMRL »SUMI M ) 
60 TO 110 
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Subroutine sumpro <lo*i*nc*max»facrl*facim*sumrl»sumim) 

DIMENSION NC(25)»FACRL<25) ,FACIM(25) 
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NC(1 )=N0 
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PDRL2»PDRL*CCRL-PDIM*FACIM(K) 
PDIM2-PDRL*FACIM(X)+P0IM*CCRL 
PDRL=PDRL2 
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10 CONTINUE 

IF (KtOUT) 250f250tl60 
15 K0UT=1 
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IF (RAS(J)) 6»3,6 
KXJ«KX-J 
SARL(XXJ)=0.0 
SAIM(KXJ)=0,0 
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SARDIF=0*01*A6SF (RAS(J)) 
bARERR*ABSF ( SARC ( J ) -RAS( J ) ) 
IF (SARDIF-SARERR) 30,40,40 
30 KARM=1 
40 CONTINUE 
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54 IER»1 
52 NP3=N+3 
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S»EXPF(S/FPN1) 
eO 208 J«3»NP3 
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356 IF(AVHB1-SK)450. 354*354 
354 B(NP3)»H(NP3)-0*B(NP3-2) 

400 IF{B(NP3) )401»550*401 

401 AVHB2=ABSF(H(NP3) /B(NP3) ) 
403 IF(SIC-AVHB2)550*450,450 



450 BO 451 J=3,NP3 

E>( J)*H( J)4R«D( J-1) 

451 E( J)=D( J)+R*E( J-1) 
IF(D(NP3) )452. 500.452 
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ITERPR IS A PRINTOUT FORMAT CHECKING CONVERGENCE FOR THE 
BANDWIDTH SOLUTION. 
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